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1 Introduction 


While wind energy has been employed for electricity production since the 1880s, it wasn’t 
until the oil crisis of the 1970s that commercial wind energy production was pursued actively 
in the United States. Wind energy use has grown rapidly since it began to be promoted as an 
alternative to fossil fuels and was accorded sponsorship by the state of California in the 1980s 
and by the Federal Government beginning in the late 1990s. Concer ns about avian and chi¬ 
ropte ran deaths caused by wind turbines emerged in the early 1990s Howell and DiDonatd. 


199l| . with widely varying e stimates of the fatality rates, and studies were mounted to assess 


these rates as early as 1998 [Smallwood and Thelandeil . 2005 1. Aggregate U. S. mortality es¬ 
timates have been reported ranging from 20p00 to 573,000 b i rds annuply Erickson et al.l . 
2001 . 2005 . Loss et ah . 2013 . Manville . 20091 Smallwood . 2013 . Sovacooi 2012 |. High profile 


lawsuits in such places as Altamont, CA (2007), Ventura, CA (2012), Nantucket Sound, MA 
(2012), Port Clinton, OH (2014) have brought the issue to national prominence. 

The naive approach to estimating turbine-related avian and chiropteran mortality— sur¬ 
veying periodically for bird and bat carcasses in designated areas near turbines at prescribed 
time intervals, and scaling the counts by time interval and study area— leads to grossly 
distorted estimates, for a variety of reasons. Some carcasses will be removed by scavengers 
before the survey, for example; some carcasses may be present but undetected at the time of 
the survey; some fatally injured birds or bats may survive long enough to alight outside the 
study area; and carcasses may be discovered whose death arose from other causes or during 
other time periods. 

A number of investigators have developed modeling approaches leading to proposed ad- 
justment formulas in t ended to overcome the distortions a r id biase s of the naive app roach 
Erickson et ah . 1998 . Johnson et ahl . 2003, Shoenfel^ . 2004J . Pollockl . 2007, Hnsol . 2011 |. each 


*Robert L. Wolpert (wolpert@stat.duke.edu, +1-919-684-3275) is Professor of Statistical Science 
and Professor of Environmental Science and Policy, Nicholas School of the Environment at Duke University, 
Durham NC, USA. 
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embodying slightly different assumptions about the processes affecting carcass discovery. The 
wide variability of these estimation formulas leaves practitioners uncertain which of them (if 
any) to use. Here we explain the assumptions that underlie four commonly used estimation 
formulas, illustrate when each is appropriate and how they differ, and propose a new model- 
based Avian and Chiropteran Mortality Estimator called “ACME” that extends all four of 
them and introduces three new features to improve the reliability of mortality estimates: the 
diminishment of Field Technician (FT) discovery prohciency as carcasses age; the reduced 
rate of scavenger removal as carcasses age; and the possibility that some but not all carcasses 
present but undiscovered by FTs in one search may be discovered in a later search. 


2 The Model Underlying the New Estimator 


Suppose that carcasses arrive in a Poisson stream with intensity m{t) that varies slowly with 
time t and that they are removed (principally by scavengers) independently after random 
times Tj with complimentary CDF F{t) = P[rj > t]. Suppose too that held technicians 
(FTs) mount blinded searches at a sequence of times T, at similar intervals A = [Tj — Tj_i], 
and that the probability that a carcass of age r will be discovered by an FT in such a search 
is S'(r) (which may depend on the carcass age r, but we are assuming for now that discovery 
is statistically independent of the scavenging removal process). Let Cj denote the (random) 
number of carcasses actually discovered in the search at time Tj. Then the expected number 
of carcasses that arrive during the period and are discovered at time Ti is 

c° := f ' m{t)F{Ti - t)S{Ti - t) dt. 

Some existing mortality estimators (see Sections fl2.lll2.2p i embody the assumption that all 
carcasses that arrived prior to the previous search at time Tj_i will have been removed by 
scavengers or discovered and removed by an FT in that earlier search, leaving none to “bleed 
through” from earlier periods to be removed or discovered in the current search at time Tj. 
Under that assumption, c° would be the expected count E[Uj]. Other mortality estimators are 
based on a different assumption— that undiscovered and unremoved carcasses from earlier 
periods remain discoverable, so that Ci may include both “new” carcasses from the current 
period and “old” ones that arrived during earlier periods. For k > 1 the expected number 
discovered at time Ti that arrived during the kth previous period but were undiscovered in 
k previous searches would be 



m{t)F{T,-t)S{T,-t) J] [1-^(T,_„ 


t)] dt 


and the total expected carcass count for the ith search would be E[Uj] = Ci := J2k>o 

Evidence (see Section ([5])) suggests that both the assumption that all carcasses bleed 
through for later discovery, and the assumption that none do, are wrong. We here introduce 
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an intermediate possibility: that some fraction 0 < IB < 1 do bleed through at each search, 
leading to expected carcass count 


oo rTi_k 

E[a] = Q;=5^53M m{t)F{T,-t)S{T,-t) J] [l - S{T,_^ - t)] dt. (1) 

k=0 0<n<k 

For slowly-varying m{t) ^ m, this leads to a maximum likelihood estimate for the mean 
total mortality m* = ^ m{t)dt ~ m/j in period (Tj_i,Tj] of 

M* := Q/Rt, (2a) 

the carcass count Cj inflated by a factor 1/R* given by the inverse of the “reduction factor” 


F(T.-t)S(T,-t) J] [1 - S(r._„ - «)] * 


(2b) 


k=0 




0<n<k 


(so-called because on average the count Cj ~ MiR* will be the mortality Mj reduced by the 
factor R*). For similar search intervals R ~ I, the fcth term in this sum for k > 1 represents 
carcasses that arrived between kl days and {k + l)I days before the end of this search period, 
were unremoved by scavengers over that entire period, were undiscovered and yet remained 
discoverable in k consecutive searches, and were finally discovered at time T,. This will be a 
rare event unless kl is quite small, so only a few terms of this sum are typically sufficient to 
achieve accuracy within a few percent. Simple approximations and truncation error bounds 
for them are given i n Section fl3.ip . 


Shoenfeldl 2004l | describes as periodic those estimators (including his own) based on 


the premise that all the undiscovered and unremoved carcasses remain discoverable, and 
the assumption that consecutive periods are similar. Our proposed estimator, intermediate 
between the periodic ones that assume 100% bleed-through and the aperiodic ones that 
assume 0%, might be described as partially-periodic. 


2.1 Special Cases Sz. Previous Estimators 

Before turning to the general case, consider first the simple situation with constant removal 
rate (or hazard) [—F'/F]{t) = r and constant search prohciency S'(r) = s. Under this 
assumption that the scavenger removal rate and FT discover probabilities do not depend 
on carcass age r, the removal times must follow the exponential distribution r r\j Ex(r) 
with survival function F{t) := P[r > t] = exp(—rf) for f > 0 and mean removal time 
t ;= E[r] = 1/r. In that case, for constant inter-search intervals R = /, the reduction factor 
fl2b|) simplihes to a geometric series. 
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1 Ai—k)I 

R* := 7 23^ / exp ( - r{il - t))s{l - s)^ dt 

7,_Q ^ {i—k—l)I 

s - 1] 


k=0 -Jii-k-l)! 

.?■/ 11 ^ ^ j-gj/t _ 


r/[e"^ - ®(1 - s)] l[eRi - ^{1 - s)] 

In the case of zero bleed-through, = 0 and ([3]) leads to the estimator 


( 3 ) 


M[ = 


I a 


s t [1 — e~R^] ’ 


(4a) 


th at introduce d bv iPollockI j2007l | (under exponentially-distributed persistence). 

Husol |201l| introduced a similar estimator that differs in replacing the term [1— 


by min (0.99, [1 — e The two are identical whenever (as usual) search intervals / are 

shorter than the mean removal times i times a factor of log 100 ~ 4.6, for then [1 — < 

0.99 (otherwise Huso’s estimator is up to 1% higher than Pollock’s Mf). 

In the cas e of fu ll bleed-through, IB = 1 and ([3]) gives the “periodic” estimator introduced 
by Shoenfeld [2004 1, 


Mf = 


iq 

s i 


J/t — 1 

g//t _ ]_ 


(4b) 


Finally, setting fB = 1/(1 — s) gives 


MF = 


iq 

s i 


the steady-state estimator introduced bv lErickson et ahl [1998 


(4c) 


2.2 Comparing Current Estimators 

All four of the estimators Mf, Mf, Mf, and special case^of ([3]), for specihc values 

of 03. Always 


[1 - = Mf <Mf <Mf <M. 


H 


(5) 


^For u nusually long search intervals h > 4.6t then is up to 1% higher than special case of 
M*. Also PollodJ ’s estimator is not limited to exponentially-distribu ted removal tim es r with constant 


removal rate r = 1ft, although the method commonly used to estimate t [Erickson et all l2008l §2.6 & §3.2] 
is the MLE for that case and is badly biased for heavier-tailed distributions. 
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so all four estimators are within 5% if J > 3t and within 58% for I > t. Under the 
assumptions of constant removal rate —F'/F = r and constant searcher prohciency S = s, 
the proposed new estimator M* of ([2]) also lies in the interval [Mf', ^1^] 0 < 03 < 1. 

Differences among the estimators will be substantial for shorter search intervals, however. 
For example, for search intervals substantially shorter than the mean scavenger removal time, 
/ -C t and so 


Mf > Mf > {i/I)Mf » Mf, 

and it will be important to assess bleed-through rate 03 accurately. And, if the assumptions 
of constant removal rates and search prohciencies are incorrect, then the estimators may 
agree with each other but all be badly biased. 


3 Variable Search Proficiency and Removal Rates 

Both the assumptions of constant removal rate and of constant search prohciency, irrespective 
of carcass age, appear inconsistent with the observations presented in Section ([5]). In this 
section we show how to go beyond those assumptions. 

3.1 Diminishing Proficiency 

For many data sets the search prohciency S{t) appears to diminish with increasing carcass 
age t. In Section ([5]) it is shown that the data are ht well by an exponentially decreasing 
success rate 


S(t) = exp — a — bt) (6) 

for parameters a, 6 > 0 (logistic models gave very similar results). With this modeling 
choice, and for equal search intervals A = / (say, with searches at times Tj = iJ), the ACME 
estimator M* and reduction factor R* of ([2]) take the form 

M* := a/i?*, (7a) 


with 


R- 



—a—xbl _ ^—a—{x—n)bl 

0<n<k 


dx 


OO 




(7b) 


whose kth term represents the fraction of carcasses that arrived in the search period 
ending at that are discovered at time R. Of particular importance (see Section (j4])) is 
the hrst of these 
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(7c) 


T* = f F{xl) e-"'’' dx, 

Jo 

the fraction of carcasses discovered at the search ending the interval in which they arrived. 
Each is expressible as the sum of 2^ terms of the form 


;= [' F{{k + (8) 

Jo 

for suitable nonnegative integers m,n that can be enumerated recursively: beginning with 
{k,m,n) = (0,1,0), each entry {k,m,n) generates at the next level (A: + l,m,?7,+ 1) and 
{k + l,m + l,n + k + 1). The hrst few terms are 


T* 

-'o 

T* 

T* 

^3 


Qoio 

Qlll + Qi21 

Q212 T Q223 Q222 T Q233 

Q313 + Q325 + ^324 + ^336 + ^323 + ^335 


+ ^334 + Q 


★ 

346 


(9) 


The truncation error from using only the first N terms 0 < /c < iV of the inhnite sum in 
fiTbll is bounded by 


0 < Truncation Error < F{N I) 

For the examples presented in Section ([5]) , the truncation error bound is about 1% of R* 
with N = 3 terms, and about 0.1% with N = 5 terms. 


^—a—N bl 


A (1 -e-“)^ 


&/ V (1 - 


( 10 ) 


3.2 Persistence Distributions 


Bispo et ah 2013alf9 | found (and we verify in Section fl5.ll) below) that log normal, log logis¬ 


tic, and Weibull distributions with decreasing hazard functions all fit empirical persistence 
data quite well, and that exponential distributions did not. Here we take the Weibull distri¬ 
bution, parametrized in the form 


F{t) := P[r > t] = exp ( — (pt)“), f > 0 (11) 

for rate p > 0 (in units of day~^) and unitless shape parameter a > 0. For this distribution 
the key quantities Q\rnn froni ([8]) needed to compute B* are 
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Qkmn ■= ?B^(— f exp ( — {p{k + x)/)" —m{a + blx)— nbl) dx, 

Jo 

easily evaluated numerically using Simpson’s quadrature rule or, for the particular values of 
a = I and a = 1, available explicitly in closed form: 


Qlrnn =-IT-~ — n)bI + p/Amb) 


pi 


{a = i) 


X 


P 


2mb 


^—mb{y/kTp+p/2mbp —mb{’^J(k+l)Ip+p/2mbp 


Tip 


^[^/2mb/p{^/lM+^)) - ^U2mb/pyp{k+l)I + 


Q kmn 


mb 

^ Q—m{a+bI) 


2mb' 


p 


2mb' 


[l_g-(p+m6p] 


{p + mb)I 

where *h(^) denotes the CDF for the standard No(0,1) normal distribution. 


(a = 1) 


4 Mortality Estimates 


Point estimates like M* of fl2al) and fl7a|) a re more inforruative w hen accompanied by some 
measure of their uncertainty. For example, lErickson et ahl 1998l | recommend reporting 50% 
and 90% interval estimates for mortality. 


4.1 Interval Estimates for Mean Mortality 

In this section we will hnd interval estimates for the mean daily mortality rate ruj based on 
observed carcass counts Ci. Such an estimate is given by a pair of functions lo(c) and hi(c) 
with the property that 


P[mi e [lo(Ci), hi(Ci)] ] > 7 

for specified 7 (such as 0.5 or 0.9, per Erickson et al.l |l998 |b The common symmetric choice 
is to arrange that P[mj < lo(Ci)] and P[mj > hi(Ci)] are each below (1— 7 )/ 2 . Frequently 
in practice however mortality is low en ough (or removal is rapid enough) that observed 
counts as low as zero or one are common |Huso et ahl. 20141, motivating interest in one-sided 


interval estimates with lo(c) = 0 and P[0 < m* < hi(Ci)] > 7 . A third option is to find the 
shortest interval that captures m with probability at least 7 . 

Under the model introduced in Sections (j2l[3j) the mortality Mi in the ith search period 
(Ti_i,Ti] has a Poisson distribution whose mean is the product mj/j of the average daily 
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mortality in that period m, and the search period length li = (Tj — Tj_i). If these rates and 
lengths are nearly constant (say, rrii ^ m and It ^ I) over the period during which all the 
carcasses found at time Ti arrived, and if the model parameters determining the reduction 
factor R* of Eqn fl7bp are nearly constant, then the conditional (given m) distribution of Ci 

is 


Cj I m ~ Po 

With conjugate Gamma prior distribution m ~ Ga(^, A) (more on this below), the marginal 
distribution of carcass counts is negative binomial 


Q~NB(e,A/(A + </)) 


( 12 ) 


and the posterior distribution for m given Ci is again Gamma but with new parameters: 


m 


Ci ~ Ga(^ + Ci, A + -Rj /). 


(13) 


The Objective Bayes reference prior distribution jBerger et al.l . l2009l | for m, expressing no 


available prior or extrinsic information about it, is the improper m ~ m~^, the limiting case 
of the Gamma distribution with ^ ^ and A = 0. An alternative to Objective Bayes is to 

follow an Empirical Bayes approach |Robbinsl . 1955 . Gasellal. 1985| using the evidence about 

m reflected by previous observations of {Gj} ~ NB(,^, A/(A + R*I)) (typically this leads to 
shorter intervals, since they reflect more evidence about the average mortality rate m). It 
proceeds by making (often Maximum Likelihood) estimates ^ and A of the parameters, and 
basing interval estimates for m on these. 

The resulting posterior 7 = 50% or 7 = 90% Gredible Interval estimates for m are of the 
form [lo(Ci), hi(Ci)] with the functions lo(c) and hi(c) given by one of: 

lo(c) = 0 

hi(c) = qgaiiima( 7 , c + A, R*i) 
lo(c) = qgaiiima((l - 7)/2, c + A 


Objective Bayes, One-Sided: 
Objective Bayes, Symmetric: 
Empirical Bayes, One-Sided: 
Empirical Bayes, Symmetric: 


c -h 


R*I) 

R*I) 


A + R*I) 


hi(c) = qgaiiiina((l-h 7 )/ 2 , 
lo(c) = 0 

hi(c) = qganima( 7 , ^ c, A R*l) 
lo(c) = qgainma((l - 7 )/ 2 , f+c, _ 
hi(c) = qganima((l-h 7 )/ 2 , c, A-hR*l) 

where qgainma(x, a, b) |R Gore Teaml . 2015j | denotes the quantile function (inverse GDF) for 
the Gamma distribution. If the mortality rate m{t) varies slowly enough that it may be 
considered constant over a longer period of time including some n > 2 search intervals of 
total length J_|_ := (Tj — Ti_n), then the total number of carcasses G+ := ^Ci found in 
the n searches will again have a Poisson conditional distribution G+ | m ~ Po(R*/_|_m) and 
a Negative Binomial marginal distribution G+ ~ NB(^, A/(A -|- 7?*/+)), and the posterior 
for m will again be Gamma, m \ C+ r\j Ga(e + G+,A + R*J+). Quantiles of this Gamma 
























distribution will determine Credible Intervals for m that will be narrower by approximately 
a factor of ^/n than those of fll3l) . and so will specify m to higher precision. The assumption 
of near-constancy of m and the model parameters determining R* would be violated for 
periods long enough to include changes in season, vegetation, or migratory patterns. 

4.2 Interval Estimates for Mortality Mj 

In this section we hnd interval estimates for the number Mj of carcasses that arrived in the 
interval (rj_i,rj] based on the observed carcass count C*. These will be wider than the 
intervals for of Section fl4.1l) because the aleatoric uncertainty and variability of mortality 
events typically exceeds the epistemic uncertainty about parameter values. 

In general the Q carcasses discovered in the search at time T, may include both some 
of the Mi carcasses that arrived during the period as well as some of those that arrived in 
earlier periods. Thus there is no way of making meaningful interval estimates about Mj from 
Ci alone, without making some assumptions about either the {Mj} for j < i, i.e., about 
mortality in the recent past, or about the absence of bleed-through. 

4.2.1 Classical Confidence Intervals (fB = 0 only) 

If, despite the evidence in Section (E]) , one assumes that no carcasses from earlier periods 
are ever discovered, i.e., if 03 = 0, then Cj ~ Bi(Mj,i?*) and classical Conhdence Interval 
estimates are available for this binomial model without concern for mortality in earlier peri¬ 
ods. For example, a 90% one-sided classical conhdence interval for Mj would be [Ci,hi(Ci)], 
where 


hi(c) =inf{M > c : pbinom(c, M, R*) < 0.10} 


where pbinom(x, n. 


p) |R Core Team . 2015 


denotes the CDF for the Binomial distribution. 


4.2.2 Objective Bayes Credible Intervals (any fB) 

No simple classical conhdence intervals for Mj are available for the more realistic situation 
of 03 > 0. Again, however. Objective Bayes and Empirical Bayes credible intervals may 
be constructed for Mj based on the model of Sections ([21E])- Both Objective and Empirical 
Bayes posterior distribution for Mj, given Cj, are derived in Appendix IA.2I and presented as 


P[Mj = M I Cj = C] = c X 2 Fii-C, -M; ^-C-M;-z) (14) 

with ^ ^ for Objective Bayes or for Empirical Bayes, for specihed quantities c and 2 ; 

given in Eqns fl21bll21ap . respectively, as explicit functions of R* and Tn fro m Eqns fjTbllTB 
(here 2 ^ 1 ( 0 ., b; c; z) denotes Gauss’ hypergeometric function |NIST DLMFI. §15]). Setting 
p(m|c) := P[M = m I C = c] from ffTTj) . credible intervals for M are 
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[lo(C,),hi(Ci)l 


(15) 


with 


One-sided: 


Symmetric: 


lo(c) =0 
hi(c) = min{M : 
lo(c) = max{M : 
hi(c) = min{M : 


J2m<MM^) > 7 } 
E,„<mP(“|c) < (1 -7)/2} 
Em<MP(“|c) > (1 + 7)/2} 


while Highest Posteri or Density or HPD in tervals [the shortest possible intervals with cover¬ 


age probability 7 , see iGelman et all l2009l. §2.3] for M upon observing Ci = c are available 


by sorting the values {p(ni|c) : m > 0 } in decreasing order and identifying the smallest collec¬ 
tion whose sum exceeds 7 . Some of these distributions and intervals are shown in Figure ([3]). 
Similar Empirical Bayes results are available from Eqns fll4ll21ap with estimated hyperpa¬ 
rameters A. 

In the absence of bleed-through {i.e., 03 = 0) all found carcasses are “new” so necessarily 
Mi > Ci- It is shown in Section flA.2.2p that the number (Mj — Ci) of undiscovered carcasses 
then has the Negative Binomial distribution (Mj — Ci) | Q ~ NB(^ -|-C, (A-|-i?*/)/(A, 


so 


P[M. = M I Cj = C] = + A//)-s-" (16) 

from which credible intervals for Mi are available. For example, the one-sided Objective 
Bayes interval is [lo(Ci),hi(Ci)] with 


lo(c) = c hi(c) = c-F qnbinom( 7 , c-|- 1 / 2 , R*) (17a) 

where qninom(p, alpha, prob) |R, Core Teaml . 2015 ] denotes the quantile function for the 
negative binomial distribution. HPD regions are available with a search. 

A more direct and less model-dependent Bayesian approach to Ending the conditional 
distribution of M given C would be to begin with an improper uniform prior distribution 
for M on the nonnegative integers {0,1,... }. The posterior distribution of the unobserved 
carcass count (M — (7), after observing C ~ Bi(M, i?*), then has the negative binomial 
distribution (M* — Ci) \ Ci ~ NB(C'i -|- 1, i?*), leading to very similar one-sided intervals with 


lo(c) = c hi(c) = c-f qnbinom( 7 , c-|-1, R*). (17b) 


5 Results from Altamont 


Warren-Hicks etEl |2012l | report on data taken from January 7 to April 30 of 2011 in the 


Altamont Pass Wind Resource Area in a study of the removal and discovery rates of aging 
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bird and bat carcasses. One hundred and ten bird c arcasses ('predominantly b rown-headed 
cowbirds, Molothrus Ater, with AOU code BHCO Pvle and DeSant^ . 2014l |i and 78 bat 
carcasses of disparate species were placed by Project Field Managers (PFMs), who then 
checked every few days to conhrm whether or not each carcass remained in place. Field 
Technicians (FTs) would search for carcasses at approximately one week intervals, noting the 
species and location of those they discovered but not disturbing or removing them. Successive 
searches were conducted by different FTs who w ere unaware of any earlier c arcass discoveries. 
This “integrated detection trial” or IDT design |Warren-Hicks et ah . 20121. Chap. 2] afforded 
the possibility of exploring how removal rates and discovery probabilities may change over 
time. 


5.1 Removal by scavengers 

Figure ([T]) illustrates the removal of brown-headed cowbird carcasses by scavengers. Removals 


Carcass Removal 


03 

E 

0 

cc 


03 
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Figure 1: Empirical survival function showing removal of brown-headed cowbirds by scav¬ 
engers (solid black stair-step lines), along with best Weibull distribution £t (solid blue line) 
and best Exponential distribution £t (dashed red line). Note Weibull fits well while Expo¬ 
nential does not. 


are interval censored: we only observe the times of the last recorded discovery of a carcass’s 
presence and the first of its absence. Thus the empirical survival function in Figure ([1]) con¬ 
sists of two black stair-step curves based on the earliest and latest possible times of removal 
consistent with the observations. The best Weibull distribution fit (see Section flA.l.lD for 
derivation of likelihood function flT8|) and MLEs), 

P[r > t] = exp ( — (pt)“), a = 0.4695, p = 0.0809 day“^ 

is illustrated with the solid blue curve. Its mean of E[r] = 27.97day is nearly twice that 
(16.49 day) of the best exponential distribution £t, shown as a dashed red line. The expo¬ 
nential distribution model underestimates early removal rates and overestimates later ones. 
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The estimated shape parameter a = 0.4695 is 9.8 standard errors away from the value a = 1 
for the exponential distribution, making the exponential distribution and its assumption of 
constant removal rates entirely untenable. Best hts with log normal and log logistic were 
nearly indistinguishable from Weibull, so we present only Weibull results here. 


5.2 Search Proficiency and Bleed-through 

Figure ([2]) illustrates the search process by FTs. Short vertical dashes at the top and bottom 


Search Proficiency for: Brown-Headed Cowbird 


CO 

CD 
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Figure 2: Empirical plot of search proficiency (dashed black line) for brown-headed cowbirds, 
along with best £t exponentially-decreasing curve (solid blue line) and best £t of a constant 
prohciency (dashed red line). Search successes (and failures) are shown as whiskers at the 
top (and bottom) of the plot, respectively. Note diminishing proficiency model fits data well 
while constant prohciency does not. 


of the plot indicate the times of successful and unsuccessful searches, respectively. Dashed 
black curve indicates a nonparametric estimator of time-dependent search prohciency, a 
moving-average double-exponential window estimator with width of 5 day. Prohciency ex¬ 
ceeds 30% initially, but falls oh at about 7%day^\ 

Solid blue line shows best exponentially-decreasing ht, based on MLEs a = 0.4695, 
p = 0.0808, and = 0.9573 found by minimizing the negative log likelihood of Eqn ([6]) (see 
Section flA.1.2p b Dotted red line shows best constant-prohciency ht. 

The deviance between the proposed model and the constant-p rohciency mo del, a sub¬ 


model with 5 = 0 and = 1, is D = 22.63. By Wilks’ theorem |Wilksl. 1938 ] this would 


have approximately a distribution with two degrees of freedom if the constant-rate model 
were correct, evidently an entirely untenable supposition with P-value about 10“®. 

Carcasses were later discovered after an initial miss 9 times in this study, and after some 
earlier miss 12 times, conhrming that some bleed-through occurred. Estimated bleed-through 
rate is 03 = 95.73%. Evidence against full bleed-through 03 = 1 is not strong enough to 
reject that possibility. 
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5.3 Mortality Estimation at Altamont 

With the parameter estimates 


d = 0.4695 p = 0.0809 day^^ a = 1.0322 b = 0.0706 day~^ % = 0.9573 

for the Weibull removal distribution [a, p), exponentially falling search prohciency (a, b), and 
bleed-through rate (iB) (see Section flA.l|) i. we can use fl7c|) and a hve-term approximation 
to fl7bD to evaluate the Reduction Factor R* for future searches at 7 day intervals and the 
fraction of “new” carcasses Tg found in each search: 


R* = 0.2496 Tg* = 0.1740. 

This suggests that about a quarter of the carcasses are discovered eventually, 17% in the hrst 
search after arrival and the rest following bleed-through. This leads to the ACME adjusted 
mortality estimate 


M* = Ci/R* = 4.01 X Ci 


for a seven-day interval ending in a search at which Ci brown-headed cowbird carcasses are 
discovered. From the same data and parameter estimates we can End reduction factors for 
other possible search interval lengths. For example, R* = 0.14 and M* = 6.9 x Ci for I = 14- 
day searches, while R* = 0.47 and M* = 2.1 x Ci for I = 2-day searches and R* = 0.57, 
M* = 1.8 X Ci for daily searches. 

Figure ([3]) shows Objective Bayes posterior distributions (see Eqns ffmi21bj) l for the 
Brown Cowbird mortality Mj at Altamont in a 7-day search period in which Cj carcasses 
were discovered for a few small values of C*. Also given in the hgure legends are point esti¬ 
mates M* = Ci/R*, posterior means and 50% and 90% Objective Bayes posterior HPD 
interval estimates Ro and /go, derived in Section flA.2j) . These are also indicated in the hgure 
by vertical arrows at M* and and by large red squares and hlled blue disks illustrating 
/50 and /go, respectively. 


6 Discussion 

Commonly-used existing estimators give similar results if search intervals R are much longer 
than the typical time ti carcasses remain unremoved by scavengers, but differ drastically for 
more frequent searches because some of these estimators assume that undiscovered carcasses 
may remain from one search period to the next and some do not. Even when they agree they 
may be biased by disregarding the diminishing removal rate (by scavengers) and discovery 
prohciency (by Field Technicians) as carcasses age. 


13 










2 3 4 5 

Mortality m 


■ 

■ ■ 


C| = 1 

A 

■ 

Mi = 4.01 
Mi=5.66 


■ 

■ l50 = [1.5] 


• l9o-[0,12] 

• 

• 

N 


• 



• 

• 

Mi • 



• 

• 


6 8 
Mortality m 


10 


12 


(b) 




Ci = 2 

■ ■ 

■ 

Mi = 8.01 

• 


Mi = 9.67 

■ ■ l50 = [4,10] 

• 

A 

• • l90=[1J8] 

• 

• 

• 

Mi 

• 



• 


Mi • 

• 


• 

o \ 


0 


10 

Mortality m 


15 


20 


(c) 


Ci = 5 

. 1^ = 20.03 

g Mi = 21.69 
• ■ Iso = [13, 24] 
•• l9o = [8,35] 


Mi 


Mi 


W 


10 


20 

Mortality m 


30 


40 


(d) 


Figure 3: Objective Bayes posterior distribution of mortality Mi for brown-headed cowbirds 
using 7-day search intervals for carcass counts Ci = 0,1,2,5 in panels (a), (b), (c), (d), 
respectively, based on Eqns ffT^I21bj) . Large red squares show 50% HPD credible intervals, 
hlled blue disks show 90% intervals. Downward arrows indicate ACME estimates M* = 
Ci/R* and Objective Bayes posterior means M^. 
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This work presents a new estimator called ACME (an acronym for Avian and Chiropteran 
Mortality Estimator) that includes many existing estimators as special cases, but that ex¬ 
tends them in three ways: it reflects diminishing removal rates; it reflects decreasing discovery 
proficiency; and it allows for an arbitrary rate of “bleed-through” of carcasses that arrived 
before the current search period began. It also includes interval (as well as point) mortality 
estimates. 

Mathematical formulas and computational methods are derived and presented here for 
both the initial problem of estimating the model’s five parameters on the basis of field 
discovery trials, and the continuing problem of constructing point and interval estimates for 
mortality on the basis of these parameter estimates and subsequent observed carcass counts. 


Data Accessibility 


A software package acme in the open-source R computer environment |R Core Teaml . 12015 


is available at GRAN for finding maximum likelihood estimates of the model parameters and 
for evaluating the ACME estimator M*, to make use of this estimator more accessible. Data 
used in preparation for this paper are included in that package. A guide to the design of 
integrated discovery trials suitable for supporting inference about the diminishing rates of 
discovery and removal (often unavailable from current discovery trial protocols) is also under 
development. 
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A Appendix: Computational Details 

A.l Parameter Estimates 

In this section we construct maximum likelihood estimates from Integrated Detection Trial 
(IDT) data for the five parameters (a, p, a, b, 03) needed for the model of Sections (I2l[3]) to 
support point estimates M* := Ci/R* of Eqn ([^ and interval estimates [lo(Ci), hi(Ci)] of 
Section fl4.2l) for mortality Mi. 
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A. 1.1 Removal 


Persistence times in this model have the Weibull distribution ffTTl) with P[r > t] = exp ( — 
for t > 0, depending on the two parameters a and p. Carcass placement times to are 
known, but removal times tr (by scavengers) are generally not observed. The data available 
from an IDT bearing on {a, p) from the /cth carcass consist of its placement time fg, the last 
time fp > of its known presence from discovery by either a FT or PFM, and the hrst time 
> tp of its conhrmed absence by a PFM (or = cxo if it remains present throughout the 
trial). The negative log likelihood function on the basis of these interval-censored data is 


4em(a, P) = -Y1 } 

k 

= ^ log |i - 

k k 

The MLEs presented in Section 05.1 p are the minimizing values {&,p), easily found by a 
numeric search, along with approximate standard errors from the inverse Hessian. 



A. 1.2 Discovery 

The probability of discovery of a f-day-old carcass present at an FT’s search is given in ([H]) 
as S(t) = exp (^ — a — bt), depending on the two parameters (a, b). 

Again denote by to the placement time for a particular carcass (say, the kth) and by tp 
the last time it is known to be present. Let mo := minjn : > to} and m* := max{n > 

'■ Tn < tp index the hrst and last FT searches at which the carcass is present, and let 
m* := max{?7, > mg : = 1} index the last successful search (or m* = mo if it is never 

discovered). Introduce the short-hand notation pn{a,b) := exp ( — a — b(Tn — to)) for the 
probability of discovery at the nth search, for mo < n < m*. For a carcass that arrived in 
an earlier search period to be discovered now it must have been undiscovered and also “bled 
through” at each previous search. Set Dn = 1 for a successful discovery and Dn = 0 for a 
failure. Then the probability of the observed sequence of successes and failures for the kth 
carcass, as a function of (a, 6, 03), is the sum over all possible indices m of the last search 
time Tm at which the carcass bleeds through. 


CiJa,b,<S) = {l-<B) n Pn{a,bp~{l-p„{a,b)y-‘’" 

mo<n<m 

p4a,b)^-{l -p„(a,6))i-^". 

mQ<n<m* 

The negative log likelihood contribution for all carcass combined is the sum 


(19a) 


£disc(a,6,®) = 5^-log£^i,,(a,6,®). (19b) 

k 
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The MLEs presented in Section fl5.2p are the minimizing values (a, b, t8). 


A.2 Posterior Distribution of Mortality 


In this section we consider the posterior distribution of the mortality Mj in a fixed period 
(Tj-i, Ti] of length li = (Tj —Tj_i) days, conditional upon the observed count Ci in the search 
at time Tj, in order to hnd interval estimates for Mj. To make the notation less cumbersome 
we omit the subscripts “i”. 

The total number C of carcasses discovered in the search will in general be a sum C = 
C'new + C'oid of “new” carcasses that arrived during the current interval and “old” ones that 
arrived in earlier periods, but were undiscovered and unremoved in earlier periods. In this 
model the mortality M ~ Po(ml) in a particular search interval has a Poisson distribution 
with uncertain mean ml for a daily average rate m > 0 which varies sufficiently slowly 
from one interval to another that we may treat it as constant over the arrival times of all 
the carcasses discovered in a particular search. We employ a Gamma prior distribution 
m Ga(£, A) for m , usually with the Objective Bayes prior parameters ^ = 1/2, A = 0 


Berger et all. 12009 


Each of the M carcasses that arrive during the period has probability Tq of being discov¬ 
ered in the current search, probability (i?* — Tq) of being discovered in some future search, 
and probability (1 — R*) of never being discovered. Thus the model may be described: 


m ~ Ga(^, A) 

M ~ Po(m/) 

G_|M~Bi(M,T*) 

0,ew ~ Po(mTo*J) 

C = Gnew + Cold 


Average daily mortality 
From all previous periods 
Mortality this period 
New count, conditional on M 
New count, marginal 
(new -|- old), indep. 


where R* > Tq are given in Eqns (frbj lTcll. 


A.2.1 Possible bleed-through (03 > 0) 

First consider the case where R* > Tq , and in particular 03 > 0, so the mortality M may 
take any nonnegative integer value— even M < C, since some or even all of the C discovered 
carcasses may have arrived in earlier search intervals. Summing over the possible number 
X = Gne„ of new carcasses and integrating over the uncertain mean daily mortality m. 
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(C-x)! 
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'j^*^M—xjC+M- 


[X + {R* 


T* + l)/]«+c^+^- 


= c X 


= C X 


CAM 

E 


r(e + c + M - x) 
^ (C — x)!(M — x)!x! 


x=0 

r(e + c + M) 

C\M\ 


2Fi(-C',-M;l-e-C'-M; -z) 


(20a) 


where 2 Fi{a,h]C] z) is Gauss’ hypergeometric function NIST DLMFI . §15] and where 


c = 


_ T*[X + {R* - T* + 1)/] 
r(0[A + {R* - T* + 1 )IY+C+M ^ - (1-T*)(i?*-T*)/ ■ 

The induced marginal distribution of G ~ Po (mi?*J) is negative binomial, 

P|C = CJ = ‘tkIct 


(20b) 


Dividing fl20ap by fl20bp gives the conditional distribution for mortality M given a carcass 
count of C\ 


P[M = M \ C = C] = cx 2 Fi(-G, -M; 1 - ^ - G - M; -^) 
with c and z given by 

_ r(e + G + M)(A + i?*J)«+^(i?* - To*)^(l - Tq*)^ 

r(^ + C) M\ (fi*)0|A + (fl* - To* + l)/|S+c+" 
r„» [A + (g - r„- + 1)1] 

(i-r*)(fl*-r„*)/ ■ 

For the Objective Bayes reference values ^ ^ and A = 0 the distribution is again given by 

ffTT|) . but c and 2 ; are a bit simpler and don’t depend on the search interval length I: 

_ r(i + G + M){R*)^R* - T*f{l - T*)^ 

r(i + c) M< (fl* - r„* + i)i+c+« 

TsiR- - T,; +1) 


(B 


(21a) 
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This is well-defined even though, by (IT^ . the marginal predictive distribution of C is degener¬ 
ate for A = 0. Objective Bayes 1007 % credible intervals for M are presented in Section fl4.2.2p 
and illustrated in Figure (E]), based on the conditional distribution of M (for specified C) 
given in Eqns fimi21bp . 

A.2.2 No bleed-through (03 = 0) 

For the remaining case of R* = Tq where only “new” carcasses can be found, the sum of 
Eqn fl20ap reduces to the single term x = C < M, so only values M > C are possible and 
for these Eqn (lldh becomes: 


P|M = M I C = C| = + + A//)-«-", (Ha 

The number M — C oi undiscovered carcasses will have a negative binomial conditional 

distribution (M — C) | C ~ NB(,^ -|- (F, (A -|- R*I)/{X + /)) or, for the Objective Bayes case 

of A = 0, (M — C) I O ~ NB(C' -I- |, R*), justifying the interval estimate for M given 

in Eqn fll7ap . 
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